arXiv: 1506.09204v4 [cond-mat.str-el] 21 Dec 2015 


Quantum spin Hall effect on the porphyrin lattice 
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We study the physics of quantum spin Hall (QSH) effect and topological quantum phase transition 
on the porphyrin lattice. We show that in a special limit the pristine model on this lattice reduces 
to the usual topological insulator (TI) thin film model. By incorporating explicit time reversal 
symmetry (TRS) breaking interactions, we explore the topological properties of this model and 
show that it exhibits an ordinary insulator phase, a quantum anomalous Hall (QAH) phase, and a 
quantum spin Hall (QSH) phase. 
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I. INTRODUCTION 

The study of topological insulators have become a 
ubiquitous study in condensed matter physics. These 
materials have captivated considerable attention of 
researchers 1-14 . They involve materials such as Bi 2 Se 3 or 
Bi 2 Te 3 . One of the special properties of these materials 
is that they have an electronic bulk structure with a fi¬ 
nite gap separating the conduction band from the valence 
band, but their edges (for two-dimensional (2D) TIs) or 
surfaces (for 3D TIs), however, have gapless states, which 
are protected by TRS. Provided that the system remains 
time-reversal invariant, these gapless states are robust. 
The bulk Hamiltonian near the energy gap closing points 
(Dirac points) replicates a 2D Dirac Hamiltonian, which 
describes the surface states. Interesting electronic trans¬ 
port properties emerge when TRS is explicitly broken. 
There are different ways to induce an explicit TRS break¬ 
ing interactions. They can be induced by applying a per¬ 
pendicular magnetic field or by depositing a ferromagnet 
on the surface of a TI. These TRS breaking interactions 
lift the degeneracy of the Dirac points and introduce a 
finite gap on the surface of a TI. The upshot is that one 
observes interesting electronic transport, which include 
a half quantized conductivity, Majorana bound states, 
etc. 2 ’ 15-17 . 

However, the quantum spin Hall (QSH) phase 8 ’ 18 is 
obtained by preserving TRS by the inclusion of an in¬ 
trinsic spin-orbit coupling (SOC). This phase is reminis¬ 
cent of two copies of quantum anomalous Hall (QAH) 
phases, where each copy breaks TRS and the total sys¬ 
tem remains time-reversal invariant. It has been shown 
that QSH effect can persist even when TRS is explicitly 
broken 10,19 . These interesting observations have paved 
the way for numerous investigations and the fates of 
QSH effect in complicated systems. In the QSH phase, 
there are two pairs of gapless edge state modes counter- 
propagating in the vicinity of the bulk gap, which can 
be captured in the tight binding limit of the low-energy 


Hamiltonian. The importance of tight binding model in 
topological insulators was first elucidated by Bernevig, 
Hughes, and Zhang BHZ 8 on a square lattice. Follow¬ 
ing the same line of study, Huichao Li et al. 3 coined 
an alternative tight binding model for TI thin film on a 
square lattice. In recent years, topological insulators have 
also been investigated on a honeycomb lattice, 4 ’ 10 ’ 18 ’ 20 , 
kagome lattice 21 Lieb and perovskite lattices 22 , and py- 
rochlore lattice 23 . 

The purpose of the paper is to provide another lattice 
geometry in which the physics of QSH effect and topolog¬ 
ical quantum phase transition are manifested. We inves¬ 
tigate a tight binding model on the porphyrin lattice 1 , 
in which TRS breaking through the Zeeman effect re¬ 
quires the nearest neighbour sites to have a complex am¬ 
plitude, reminiscent of the Haldane model 2 ; see Fig. (1). 
By incorporating an additional phase (gauge) into the 
complex amplitudes we show that in a special limit and 
with a specific gauge choice, the pristine model reduces 
to a time-reversal invariant Hamiltonian, which in turn 
reduces to the usual TI thin film model in the contin¬ 
uum limit. We also show that this model possesses dif¬ 
ferent Chern numbers in different regimes of the param¬ 
eter space just like in real systems. We further explore 
this model by incorporating spin degrees of freedom and 
explicit TRS breaking couplings, such as an off-diagonal 
coupling (electrostatic potential), tu, a magnetic field , h, 
and a Rashba SOC, t/j. We explicitly map out different 
phase diagrams at the gap closing points. 

II. TIGHT BINDING MODEL 

Lattice regularizations of the low-energy Hamiltonian 
is an expedient way for studying the topological proper¬ 
ties and the edge states of a TI in the entire Brillouin 
zone. In this paper, we study the tight binding model 
of porphyrin thin films. 1 The tight binding model can be 
generally written as 
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Lattice Hamiltonian of this form is usually considered as 
a toy model used to understand the physics of QSH ef¬ 
fect in real materials . 3,8,22,23 Here, cj (ci^) creates (an¬ 
nihilates) an electron with spin a (a =f, 4 _) at site 17 
r = ±1 denotes sublattice A and B respectively; and 
where i = x,y,z and r_ = (t x - ir y )/ 2 repre¬ 
sents the components of sublattice (layer) pseudospins. 

= — & m i = ±<£> for to = l ± 7i i2 respectively, i.e., 
if the electron makes a right or left diagonal turn to get 
to the second bond. The coordinates 71 2 are defined 
in Fig. (1). Aim = for m = l ± 71 , 2 ; where Ai ) 2 
are complex variables. The summations over (Im) and 
(( Im)) denote nearest and next nearest neighbours on a 
porphyrin square sublattice respectively; see Fig. (1). 

In the absence of the phase <I> and the spin degrees of 
freedom, the first three terms in Eq. (1) correspond to the 
model proposed by Joel Yuen-Zhou, et al., 1 in the context 
of Frenkel excitons (quasiparticles with no net charge). 
A special property of this tight binding model is that the 
nearest neighbour hopping sites in Eq. (1) have complex 
amplitudes, which are explained in Ref. [1]. The motiva¬ 
tion for the phase factor (gauge) <f> is rather technical. It 
does not, however, change the complex nature of the hop¬ 
ping amplitudes. In particular, the total phase in a unit 
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cell vanishes, reminiscent of the Haldane model. 2 The 
reason for this phase factor will become apparent in the 
subsequent sections, where we show that different choices 
of $ lead to different surface Hamiltonians, perhaps not 
different physics. The second and the third terms repre¬ 
sent the diagonal couplings between the upper and the 
lower surface states of the thin film. In addition, we 
have introduced other couplings in order to enhance the 
topological properties of this thin film. The fourth term 
is an electrostatic potential which couples (off-diagonal) 
the upper and lower surface states; the fifth term can be 
identified as the Rashba-like SOC in the thin film, and 
the last term is a staggered sublattice potential, which 
corresponds to a magnetic field in the thin film. 

Next, it is expedient to Fourier transform into momen¬ 
tum space, i.e., ci = ^k c k e_k Z - We obtain 
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where p 1 = SHAi, p 2 = fHA 2 ; and pi = 3Ai, p 2 = 3A 2 ; 
D\ and 3 denote real and imaginary parts. The Pauli 
matrices (J l act on the real spin space, while r, act on 
the upper and the lower surfaces of the thin film, with 
spin-up on the upper surface and spin-down on the lower 
one and I T is a 2 x 2 identity matrix. 

III. QUANTUM SPIN HALL STATE 

In order the see how the physics of QSH effect is man¬ 
ifested on the porphyrin lattice, we consider a special 
limit of Eq. (3) . This limit corresponds to the first three 
terms in Eq. (3) {i.e., all the terms involving the big and 


small square brackets), and we call it the pristine model. 
The first two terms (big square brackets) in Eq. (3) yield 
different functions depending on the choice of <f > . The 
gauge choice $ = 0 amounts to the Hamiltonian studied 
in Ref. [1]. We consider an alternative gauge choice of 
the pristine Hamiltonian with $ = n/2. In this gauge, 
the Hamiltonian decouples into two TRS copies, i.e.,T : 
% -j- i —y Tii, where the time reversal operator is given by 
T = icr y )C and A is a complex conjugation. Next, we 
take the limit pi = p 2 = P and pj = —p 2 = p, which 
corresponds to Ai = A£. The upshot of these limits is 
that the Hamiltonian simplifies as: 
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FIG. 1: Color online. Porphyrin lattice with two sublattices 
A and B. Nearest neighbour sites are along the diagonals 
with coordinates ±71 and ± 72 ; where 71 = (x + y )/2 and 
72 = (—x 4- y)/ 2. Next-nearest neighbour sites are along the 
horizontal and vertical directions with coordinates ±* and ±y 
respectively. 


(i) t m = 2\t± 



(in) 0 <t m < 2 |£j_| 



(U) t m = — 2|£_l| 



(iv) t m > 2 |£j_I 



+ [tm ~ tj_ (cos k x + COS ky))T z a z . 

Expanding Eq. (4) near the F point (k = 0), we see that 
the resulting low-energy Hamiltonian is analogous to that 
of a 3D TI thin film. 3- ' As the pristine Hamiltonian, 
Eq. (4) , has a different structure from previously studied 
model on the square lattice , 3 ’ 8 we will show how topolog¬ 
ical effects are manifested in Eq. (4) just like in the previ¬ 
ously studied models. As the system possesses TRS, the 
corresponding eigenvalues are two-fold degenerate. The 
topological phase transition points are located at the mo¬ 
menta where the energy eigenvalues vanish. There are 
four semi-metallic states at the time reversal invariant 
momenta (TRIM) points: Qo = (0,0) and Qi = (ir,ir) 
for t m = ±2t± respectively (see Fig. (3)) and Q 2 = (0, 7 r) 
and Q 3 = (7 r, 0) for t rn = 0 respectively . We focus on 
the points Qo and Qi, with nontrivial phase boundaries. 
The phase boundaries separate regions with QSH phase 
and ordinary insulator phase. The spin Chern number 
characterizing these different phases is given by 24 
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FIG. 3: Color online. The spin up energy bands of the tight 
binding model in Eq. (4) showing the metallic phases ( i ) and 
(ii), the QSH phase (in), and the insulating phase (iv). 


The Berry curvature Q/(fc) has the form: 
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where the summation is over all occupied valence bands 
in the first Brillouin zone below the bulk energy gap, and 
Vi = 9"H(k )/dhi\ i = x,y define the velocity operators. 
The Chern numbers can be computed separately. The 
appropriate formula is given by 25 
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where e aj g 7 is an antisymmetric tensor; n t are the com- 
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FIG. 4: Color online. The bulk energy spectrum and the 
edge states for the pristine model along the k y direction, with 
t± = 1 and p = 1 . (a) t m < 2|tj_ |; (6) t m = 2|tj_|; (c) 

t m > 2 |tj_|. 


ponents of a unit vector with, 



n z = (t m - t± (cos k x + cos k y ))r z . ( 10 ) 


For this model, all the six terms in Eq. (7) contribute 
equally. The numerical computation of the spin Chern 
numbers in the entire Brillouin zone yields 

fo, t m G ( oo, 2 |I_l|) and t m € (2|fJ,oo); 

C± = < undefined, t m = ±2|ij_|; 

[=Fl, tm e (—2 |*_l|, 2|tj_|) ; 

( 11 ) 

independent of p. The trivial Chern number corresponds 
to an ordinary insulator phase. The Chern number is ill- 
defined at the gap closing points (semi-metallic phase). 
The gap further reopens, which corresponds to a QSH 
phase with a nontrivial Chern number. The phase dia¬ 
gram is depicted in Fig. (2), where the phase boundary 
t m = ± 2 |£jJ separates different phases with trivial and 
non-trivial Chern numbers. 

The nontrivial Chern numbers in the QSH phase in¬ 
dicate that there are a pair of counter-propagating edge 
states in the vicinity of the bulk gap. We have solved 
for the edge states by diagonalizing the Hamiltonian in a 
cylindrical geometry, with periodic boundary conditions 
along the y direction and open conditions along the x 
direction 1 '. In Fig. (4), we show the plot of the bulk 
energy spectrum in the three phases shown in Eq. (11). 
The edge states clearly reveal the topological properties 
of this system. It is apparent that in the QSH regime 
(a), there are two pairs of edge states modes with oppo¬ 
site pseudospins (which are degenerate) that propagate 


in the vicinity of the bulk gap; these states are gapless, 
e(k y ) = 0 at k y = 0 , corresponding to the absolute value 
of the Chern number 26 . The bulk gap closes at ( b ) corre¬ 
sponding to the phase transition point. In the insulating 
regime, the gap reopens and the edge states participate in 
the bulk gap separation (c), with no counter-propagating 
edge states modes in the vicinity of the bulk gap. 


IV. TRS BREAKING INTERACTIONS 

In this section we study the topological phases of the 
full Hamiltonian in Eq. (3). Henceforth we work with the 
gauge 2 ' $ = 7t/2 and take Ai = A£, which corresponds 
to pi = p 2 = p and pi = —p 2 = p. We further take p as 
the energy unit. The full Hamiltonian in Eq. (3) does not 
possess any analytical diagonalization. We find that at 
TRIM points Qo and Qi the topological phase boundary 
is given by 



which is determined by the condition that the conduction 



FIG. 5: Color online. The bulk energy gap for the low energy 
bands of the full Hamiltonian at Qo as a function of h for fixed 
value of t± = 0.8, t« = 0.01. The red dash line is for t m = 1 
and the blue dash line is for f m = 2. There are two gapped 
regions separated by a gap closing point in the QSH regime 
tm < 2|t_i_| and the ordinary insulating regime t m > 2\t±\ of 
the pristine model. 

and valence bands touch each other as in the case of 
graphene 10,19 . Here, Ir = ta + til| t s = t m + 2st±, and 
s = =F at Qo and Qi respectively. The phase boundary 
at the other TRIM points Q 2 and Q 3 corresponds to 
t_ l = 0 in Eq. (12). Note that the topological phase 
boundary diverges at the phase boundary of the pristine 
model t m = ±2 1±. Now we consider the QSH regime 
t m < 2|fj_| and the ordinary insulating regime t m > 2|fj_| 
of the pristine model, with fixed value of f y, and several 
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FIG. 6: Color online. The bulk energy spectrum and the edge states for the full tight binding model along the k y direction. 
The parameters are: (a) h = 0.1 (b) h = h c (c) h = 0.8 with f m > 2|tjJ and (d) h = 0.1 (e) h = h c (/) h = 0.8 with t m < 2|fj_|- 
Other parameters are the same as in Fig. (5) with tR = 0.01. See text for explanation. 


values of tR. We find that the bulk gap as a function 
of h at Qo initially increases at h = 0, and gradually 
decreases as h increases from zero, and then vanishes 
at the topological phase boundary h c . The gap further 
reopens for h > h c \ see Fig. (5). This is not different 
from what is observed in other lattice geometries such 
as the honeycomb lattice, 4 ’ 10,18 ’ 20 , the kagome lattice 21 , 
the Lieb and perovskite lattices 22 , and the pyrochlore 
lattice 23 . 

We have solved for the edge states to substantiate the 
topological properties of these phases. In Fig. ( 6 ), we 
plot the edge states along the k y direction. In the ordi¬ 
nary insulating regime ( t m > 2 |fj_|) of the pristine model 
with t± = 0 . 8 ; t\\ = tn = 0.01 (a) — (c), we see that 
there is an ordinary insulating phase with no edge state 
modes at h < h c . The gap closes at the phase bound¬ 
ary h = h c which indicates a phase transition. As seen 
in Fig. ( 6 ), at h > h c the magnetic field dominates, we 
see that the gap reopens and there are a pair of counter- 
propagating edge state modes in the vicinity of the bulk 
gap, which indicates a QAH phase. We find that no 
QSH phase appears in this regime. On the other hand, 
in the QSH regime ( t m < 2|ij_|) of the pristine model 
with t_ l = 0 . 8 ; ty = tu = 0.01 ( d) — (/), there are 
two pairs of counter-propagating edge state modes in the 
vicinity of the bulk gap at h < h c , where the red and 
blue colours represent the edge state modes with opposite 
pseudospins. Again the gap closes at the phase boundary 
h = h c which indicates a phase transition. At h > h c , the 
magnetic field once again dominates and the gap reopens 
with a pair of counter-propagating edge state modes in 
the vicinity of the bulk gap, indicating a QAH phase. 


On the contrary, no ordinary insulator phase appears in 
this regime. The phase diagram in this regime is shown 
in Fig. (7). It is interesting to note that a QAH phase 
appears in each regime of the pristine model. We have 
checked that the continuum Hamiltonian also corrobo¬ 
rates our results. Thus, the porphyrin lattice provides 
another lattice geometry to study QSH effect and topo¬ 
logical quantum phase transition. 



FIG. 7: Color online. The phase diagram shown for fixed 
values of tii > 0 and t± > 0 with £ m < 2|fx|. The red dot 
represents the point where the boundary vanishes. 
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V. CONCLUSION 

The porphyrin lattice represents another lattice geom¬ 
etry where the physics of quantum spin Hall (QSH) effect 
and topological quantum phase transition are manifested. 
In this paper, we have explicitly shown how these physics 
can be constructed from a tight binding toy model. We 
showed that the porphyrin lattice tight binding model 
captures the topological quantum phase transitions ob¬ 
served in other lattice geometries. In particular, the tran¬ 
sitions between QSH phase and ordinary insulator phase, 
as well as quantum anomalous Hall (QAH) phase. We 
also showed that the complex hopping amplitudes and 


the gauge degree of freedom are the distinguishing fea¬ 
tures of the porphyrin lattice tight binding model as was 
first introduced in Ref. f 1]. 
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